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We investigate the electronic, magnetic and orbital properties of Lao.sCao.sMnOs perovskite by 
means of an ab initio electronic structure calculation within the Hartree-Fock approximation. Using 
the experimental crystal structure reported by Radaelli et al. [Phys. Rev B 55, 3015 (1997)], we 
find a charge-ordering stripe- like ground state. The periodicity of the stripes, and the insulating CE- 
type magnetic structure are in agreement with neutron x-ray and electron diffraction experiments. 
Ifowever, the detailed structure is more complex than that envisaged by simple models of charge 
and orbital order on Mn d-levels alone, and is better described as a charge-density wave of oxygen 
holes, coupled to the Mn spin/orbital order. 

PACS numbers: 75.47.Lx, 75.10.-b, 71.27.-ha, 71.45.Lr 



I 

o 
o 



> 

o 
m 
o 



I 

o 
o 



X 
J3 



Understanding the physics displayed by manganite 
oxide compounds Ri-ajD^^MnOa (R=rare-earth atom, 
D=divalent substituent) has stimulated much experimen- 
tal and theoretical work [1]. Changing the composition 
X, they show a variety of phenomena, such as ferromag- 
netic (FM), anti- ferromagnetic (AF), charge and orbital 
ordering (CO and 00) revealing that charge, spin and 
and lattice degrees of freedom are closely interrelated. 

In the present work, we concentrate on the half-doped 
case of Lai_a;Caa;Mn03 that has been extensively stud- 
ied [2-6]. For a; = 0, LaMnOa, there is a collective Jahn 
Teller (JT) distortion of the structure that gives rise to 
a splitting of the levels. As a consequence of that, the 
structure is orthorhombic, the Mn Cg-levels are singly oc- 
cupied and all Mn ions are Mn'^+ . For the other extreme 
doping [x = 1), CaMnOa is pseudo-cubic as the e^-levels 
are empty and Mn''+ does not couple to the lattice dis- 
tortion. In LaMnOa, the splitting of the Cg orbitals fa- 
vors the cooperative JT distortions and the appearance 
of 00. This fact allowed Goodenough [7] and Kugel and 
Khomskii [8] to explain the magnetic structure of this 
compound. Their explanation assumed unpolarized 0^~ 
ions and was based on the idea that magnetic ordering is 
dictated by the orientation of the orbitals involved. 

The regime of intermediate doping is complex, mag- 
netically and structurally, but even at the level of the 
local electronic structure there are several experiments 
that present contradictory results. Some of them claim 
a mixed valence picture of Mn'^+ and Mn*+ while others 
found features that do not reconcile with this image. For 
example, by doing Oxygen K-edge electron-energy-loss 
spectra, it was concluded that carriers in Lai_a;Srj;Mn03 
(0< X < 0.7) have a significant oxygen-p hole character 
[9]. Some x-ray absorption studies at the Mn K-edge 
in Ca-doped LaMnOs revealed a picture that does not 
match with a mixture of Mn'^^ and Mn^+ for interme- 
diate dopings [10]. On the contrary, Tyson et al. [11] 
performed a Mn K/j x-ray emission experiment in the 



same compound and found that their data was compat- 
ible with a Mn'^+/Mn'*+ mixing for intermediate com- 
positions. Similar contradictions have been found using 
neutron diffraction [12], photo-emission and x-ray spec- 
troscopy techniques [13]. These contradictory features 
motivate the need for a clarifying picture regarding the 
valence state of the Mn ions in manganite systems. One 
of the goals of the present work is to provide an ab initio 
quantum mechanical calculation to investigate the local 
electronic structure of the compound in order to establish 
the nature of the charge carriers. 

In this letter we report the theoretical finding of 
oxygen-stripes in Lao.sCao.sMnOs by means of an ab ini- 
tio spin-unrestricted Hartree-Fock (HF) [14] study [15]. 
Our results were obtained by means of the code CRYS- 
TAL98 [16] that uses a few localized basis functions per 
atom to solve self-consistently the HF equations. In the 
present work, the basis set for the different atoms are 
those optimized for previous calculations [17,18] and the 
crystal structure is taken from the experimental work re- 
ported by Radaelli et al [3]. So, the present study takes 
into account the actual JT distortions and the size of the 
atoms involved. It is expected that the periodic HF ap- 
proximation will correctly describe the physics involved 
in this system as it has been successful in describing other 
magnetic insulators as well as strongly correlated mate- 
rials [14]. 

The crystal structure for Lao.sCao.sMnOs can be in- 
dexed in space group P2i/m [3] with 3 inequivalent Mn 
atoms per unit cell. Two Mn ions are associated with 
JT distorted MnOg octahedra (that we will label as dl 
and d2). Although in a similar environment of local O 
distortions, these two sites are not identical. The third 
Mn belongs to an almost undistorted octahedra. The 
cell of the crystal structure has one axis (5) longer than 
the other two ones (a, c). We will call basal planes the 
Mn-0 planes perpendicular to the b axis. The magnetic 
ordering for this compound was coined CE-type [19]. It 
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consists of FM zigzag chains that are coupled AFM (see 
Fig. 1) 




FIG. 1. Top panel: Sketch of the orbital ordering ob- 
tained from ab initio HF calculations. Mn(dl)/Mn(d2) are 
the Mn in the distorted octahedra dl/d2. The "corner-Mn" 
is shown. Shading indicates the spin polarization value (see 
Table I). Small circles represent the oxygens (black ones have 
spin-polarized holes). Solid and dashed zigzag lines indicate 
oppositely spin polarized Mn species. The dotted line shows 
the unit cell used in all the calculations throughout this work. 
Bottom panel: Spin-density map of a basal plane (correspond- 
ing to the unit cell in the sketch). Continuous, dashed and 
dot-dashed isodensity curves correspond to positive, negative 
and zero values, respectively. The CE-type magnetic ordering 
and the spin-polaxized oxygen holes are clearly seen. Adja- 
cent basal planes are AFM coupled and display stacking with 
identical Mn charge-orbital ordering. 

As it was previously reported [17,20] for LaMnOs and 
CaMnOs within the HF approximation, Mulliken popula- 



tion analysis (MPA) gives charges for Mn and O ions that 
deviates substantially from the formal valence picture. In 
these works, MPA results for LaMnOa give Mn+^^ and 
O""'^''' to be compared with the formal values Mn+'^ and 
whereas in CaMnOa, Mn+^-^ and 0"^ "^ should be 
contrasted with Mn+* and 0~^. However, the MPA for 
the spin values do not show such deviation. In LaMnOs, 
the Mn spin was 1.98 [20] and in CaMnOa, ^ [17]; that 
is, very close to the formal valence spin values of 2 and |, 
respectively. In both cases, the spin for the O atoms was 
also close to the formal picture being mainly spin unpo- 
larized. The origin of the difference between charge and 
spin values is the strong hybridization between Mn and 
O bands with the subsequent screening of Mn charge by 
electrons in the surrounding O atoms. In the half doped 
case, we find that MPA also shows deviations from the 
formal valence picture with all Mn having about the same 
charge. For the oxygens, almost all of them have a formal 
charge close to 0~^'^, but there are some holes leading 
to 0~^-^, as we will detail further on. 

Regarding the magnetic properties of this system, we 
found that the ground state is correctly predicted to have 
the CE-type ordering in agreement with neutron experi- 
ments [19,3]. The FM state and the A-type AFM state 
both have energies that are above the CE-type state by 
at least 1.1 eV and 0.9 eV per unit cell, respectively [21]. 
Presumably, the FM and A-type states have similar en- 
ergies because the coupling between basal planes is weak. 

We found that the ground state for the half doped case 
is almost doubly degenerate. One of the solutions is vi- 
sualized in the spin-density map depicted in Fig. 1, that 
clearly shows that the Mn ions display 00. As an aid 
to understand the map, we draw in the same figure, a 
sketch of the corresponding orbital and crystal structure 
with shading showing the value of the spin polarization 
(S). Some oxygens have spin-polarized holes (see Table 
I). For the solution shown, holes are located on O atoms 
forming the long bonds of one of the distorted octahedra 
dl. The other (nearly degenerate) solution has oxygen 
holes instead along the long bonds of octahedra d2, with 
the corner-Mn orbital always in the direction of the O- 
hole [22]. So, within the HF approximation, not just the 
Mn ions but also the oxygens play an important role. In- 
deed, as seen in Table I, charge modulation occurs prin- 
cipally on the oxygens rather than on Mn atoms. 



TABLE I. Mulliken population values for the different 
atoms. The spin values for the Mn d-orbitals are: t2g=1.4, 

dj.2_j^2=0.1, d3^2_^2=0.5 for all of them *. 





Mn(dl) 


Mn(d2) 


corner-Mn 


0-hole t 


La 


Ca 


Charge 
Spin 


+2.1 
2.0 


+2.2 
2.0 


+2.1 
1.9 


-1.2 
-0.3 


+2.9 
0.0 


+1.9 
0.0 



* local system of reference on each Mn, with z a^xis along the 
direction of the zigzag. ^ The remaining oxygens have charge 
-1.6 and are weakly polarized. 
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The origin of the spin-polarized 0-holes can be under- 
stood in terms of the 00. For the analysis, we will con- 
sider only the atoms within a zigzag chain because the 
oxygens between chains are unpolarized (which is consis- 
tent with the AFM coupling between adjacent Mn). As 
can be seen in Fig. 1, the 00 around Mn(d2) resembles 
LaMnOa. Along the long bond of d2, we found no spin- 
polarized 0-holes which is consistent with Goodenough's 
ideas. However, around Mn(dl) the 00 is different. In 
the long bond of dl, both Mn and O orbitals lie along the 
same line allowing for a transfer of charge from the O to 
the Mn, and rendering the O with a hole [23] . Due to ex- 
change, this 0-hole has opposite spin to its two adjacent 
Mn ions. 

As the spin-polarized oxygens are always AF aligned 
with the neighboring Mn spin, the Mn-0 complex forms a 
state analogous to the Zhang-Rice singlet in the cuprates 
[24]: the corner-Mn with the 0-hole forms a low spin- 
state that gives a total S « | (see Fig. 2). This con- 
nects to the usual models of d-orbital occupancy alone, 
where a hybrid Mn-0 state stands in place of a Mn'*+ 
ion. Since we find S = 2 for the remaining Mn atoms, 
they match the conventional picture of Mn^^^ ions. So, 
although the chemistry is quite complex involving holes 
on some oxygens, due to the strong hybridization they 
are not independent degrees of freedom and the system 
should be viewed in terms of hybridized bands, possibly 
supporting the use of simplified single-band models. 




FIG. 2. Detail of a zigzag chain showing the low-spin state 
in manganites that is analogous to the Zhang-Rice singlet in 
cuprates. 

The breaking of symmetry due to the localization ol 
holes on oxygens in only one of the distorted octahedra 
appears to highlight the very small symmetry breaking in 
the experimental crystal structure, where slightly differ- 
ent octahedra dl and d2 have been reported [3] . Namely, 
in the a-c plane, the octahedra have long Mn-0 distances 
of 2.07 and 2.06 A and short Mn-0 distances of 1.93 and 
1.92 A, respectively. However, we performed a relax- 
ation by increasing the distance between the corner-Mn 



and 0-hole by 4 % and the energy went down by 0.5eV 
per unit cell, suggesting that a full relaxation is needed. 
It is worth asking if a further refinement can be per- 
formed taking into account the features suggested by our 
results. In fact, a recent refinement for the structure of 
Pro.6Cao.4Mn03 [12] produces local ordering not dissim- 
ilar to the results of our model. 

We now compare our results with electron-microscopy 
experiments (EM). As can be seen in Fig. 1, the CO and 
00 doubles the pseudo-cubic unit cell in one direction 
(vertical in Fig. 1) as reported in EM [4]. Note that the 
unit cell is further doubled in the horizontal direction 
due to the magnetic ordering although this can not be 
explicitly seen in EM. It is important to note that either 
of the doubly degenerate solutions (either holes on dl or 
holes on d2) produce identical periodicity for the diffrac- 
tion pattern. In a real sample, there will always be some 
perturbation to stabilize one of the solutions. 

In experiments, although La and Ca are on average 
randomly distributed, little is known about local arrange- 
ment within a range of a unit cell. We investigated this 
issue by checking the influence of different distributions 
for La and Ca within the planes just above and below the 
basal ones. We found that if the concentration of Ca is 
higher in the region above a particular octahedron (thus 
enhancing the symmetry breaking between dl and d2), 
then holes locate in that octahedron and the energy is 
lowered by 3.8 eV per unit cell. This result suggest that 
the local environment is important and should be taking 
into account in future refinements. 
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FIG. 3. DOS projected onto the atoms of a zigzag chain. 
The dotted line marks the Fermi energy. Indicated between 
parentheses is the octahedron to which each atom belongs. 

Our picture is supported by an analysis of the density 
of states (DOS) (Fig. 3). In the valence band, the oxy- 
gens on dl show a net down polarization (yet the other 
oxygens are almost unpolarized) and the spin-polarized 
contribution of each type of Mn is essentially up. The 
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top of the valence band is mainly populated by oxygen 
bands showing the charge-transfer insulator character of 
this system. The spin-polarized holes just at the bottom 
of the conduction band indicate that the origin of the 
gap is probably the localization of the charge produced 
by the oxygen-stripes that develop in the system, as a 
charge density wave (CDW). 

In conclusion, we performed an ab initio study for the 
half doped case of Lai_a;Ca3;Mn03 considering the ex- 
perimental structure [3,19] and working within the HF 
approximation. In contrast to the conventional model 
of CO that produces insulating behavior by the alter- 
nation of Mn'^+ and Mn'*+. our model yields an insu- 
lator because of ordering of 0-holes. Nevertheless, the 
spin character is exactly as predicted by the conven- 
tional models, because each 0-hole is bound into a lo- 
cal low spin state with the neighboring corner-Mn ion. 
The essential physics might be described with an effec- 
tive Hamiltonian involving only Mn bands, but this issue 
requires further investigation [25]. Interestingly, we find 
that this hybridized low-spin state further orders in a 
broken-symmetry CDW with AFM order. The unit cell 
and the periodicity of the oxygen-stripes are in agree- 
ment with neutron diffraction experiments. However, the 
0-holc density is substantial and not obviously compat- 
ible with the very small difference between dl and d2 
octahedra reported in Ref. [3]. Since the (La,Ca) disor- 
der will favor locally broken symmetry states, it is likely 
that these short scale fluctuations will restore the aver- 
age symmetry of dl and d2. Future work should allow 
the structure to relax within a theoretical calculation or 
to get a further refinement of the experimental crystal 
structure for Lao.sCao.sMnOa. 
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